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Abstract — Dynamic  Chest  Image  Analysis  aims  to  develop  model-based 
computer  analysis  and  visualization  methods  for  showing  focal  and  general 
abnormalities  of  lung  ventilation  and  perfusion  based  on  a  sequence  of  dig¬ 
ital  chest  fluoroscopy  frames  collected  with  the  Dynamic  Pulmonary  Imag¬ 
ing  technique  [18,5,17,6].  We  have  proposed  and  evaluated  a  multiresolu- 
tional  method  with  an  explicit  ventilation  model  based  on  pyramid  images 
for  ventilation  analysis.  We  have  further  extended  the  method  for  venti¬ 
lation  analysis  to  pulmonary  perfusion.  This  paper  focuses  on  the  clinical 
evaluation  of  our  method  for  perfusion  analysis.  Three  clinical  cases  are  in¬ 
cluded  to  illustrate  the  effects  of  contrast  media  in  perfusion  analysis  and  a 
dozen  of  clinical  cases  without  using  contrast  media  are  used  for  evaluation. 
Our  clinical  evaluation  shows  that  our  method  is  effective  for  pulmonary 
embolism,  demonstrating  consistent  correlations  with  computed  tomogra¬ 
phy  (CT)  and  nuclear  medicine  (NM)  studies.  This  performance  is  the  con¬ 
sequence  of  the  idea  that  the  cardiac  information  recorded  in  the  perfusion 
image  sequence  may  be  utilized  to  accelerate  pulmonary  perfusion  analysis 
and  improve  its  sensitivity  in  detecting  pulmonary  embolism. 

Keywords — Chest  images,  dynamic  chest  image  analysis,  pulmonary  per¬ 
fusion,  perfusion  model,  pyramid  images. 


I.  Introduction 

Lungs  take  air  in  order  to  absorb  oxygen  from  air  into  blood, 
thus  both  air  and  blood  flow  into  the  lung  are  prerequisites  to 
exchange  of  oxygen  from  air  into  blood  in  the  lung.  Inadequate 
function  of  the  lung  may  be  due  to  failure  in  ventilation  or  per¬ 
fusion  among  other  factors.  In  order  to  detect  abnormalities  in 
lung  ventilation  and  perfusion,  a  number  of  functional  and  mor¬ 
phological  studies  are  conventionally  used,  such  as,  ventilation 
isotope  scan  and  perfusion  isotope  scan.  However,  they  can  only 
provide  a  static,  coarse  geographic  2D  distribution  and  also  have 
the  disadvantage  of  using  radioactive  isotopes. 

Chest  X-ray  is  the  primary  imaging  method  for  the  diag¬ 
nosis  of  pulmonary  disorders.  Automated  analysis  of  chest 
X-ray  images  was  one  of  the  first  areas  to  receive  atten¬ 
tion  ( e.g .,  [15,2,3,1,4]).  However,  previous  work  is  mostly  re¬ 
stricted  to  a  single  chest  image  and  limited  to  using  spatial  infor¬ 
mation  for  diagnosis.  The  information  about  pulmonary  func¬ 
tion  (ventilation  and  perfusion)  that  may  be  gleaned  from  a  sin¬ 
gle  chest  X-ray  is  rather  limited,  but  it  is  evident  that  for  effec¬ 
tive  diagnosis,  the  function  of  lungs  must  be  carefully  examined. 

In  recent  years,  functional  imaging  has  become  increasingly 
prominent  as  an  important  new  frontier  in  medical  imaging  sci¬ 
ences.  Turku  University  Central  Hospital  has  taken  the  initia¬ 
tive  with  Dynamic  Pulmonary  Imaging  [18,5,17,6].  A  novel 
image  acquisition  system  has  been  developed,  which  can  grab 
a  sequence  of  digital  chest  fluoroscopy  frames  at  the  sampling 
frequency  of  25  Hz  in  a  short  period  of  time  (typically,  4.4  sec¬ 
onds).  X-rays  are  attenuated  much  less  in  air  than  in  blood  and 
soft  tissue;  consequently,  the  average  pixel  intensity  of  an  area  in 
the  lung  field  varies  over  time.  This  variation  (called  an  observa¬ 
tion)  reflects  the  air  and  blood  change  in  the  corresponding  2D 


projectional  area  of  the  lung  when  the  patient  breathes  naturally, 
that  is,  what  can  be  observed  is  a  mixture  of  the  information 
on  both  ventilation  and  perfusion  disturbed  by  noise.  This  tech¬ 
nique  opens  a  new  opportunity  for  pulmonary  ventilation  and 
perfusion  examination  without  the  use  of  radioactive  isotopes. 

Dynamic  Chest  Image  Analysis  aims  to  develop  model-based 
computer  analysis  and  visualization  methods  for  showing  fo¬ 
cal  and  general  abnormalities  of  lung  ventilation  and  perfusion 
based  on  a  sequence  of  digital  chest  fluoroscopy  frames  col¬ 
lected  with  the  Dynamic  Pulmonary  Imaging  technique  with¬ 
out  using  radioactive  isotopes.  We  have  proposed  and  evaluated 
a  multiresolutional  method  with  an  explicit  ventilation  model 
based  on  pyramid  images  for  ventilation  analysis  [10,13].  We 
have  further  extended  the  method  for  ventilation  analysis  to  pul¬ 
monary  perfusion  [9],  where  the  model  parameters  are  assigned 
with  new  medical  meanings  to  form  a  perfusion  model  and  the 
pyramid  structure  is  truncated  to  avoid  perfusion  artifacts. 

In  this  paper,  following  the  definition  of  perfusion  signals  in 
Section  II,  we  shall  review  our  model-based  perfusion  analysis 
with  truncated  pyramid  images  in  Section  III.  In  order  to  accel¬ 
erate  pulmonary  perfusion  analysis  and  improve  its  sensitivity 
in  detecting  pulmonary  embolism,  in  Section  IV,  we  present  a 
simple  yet  accurate  approach  to  extract  cardiac  systolic  and  di¬ 
astolic  phases  from  the  heart,  so  that  this  cardiac  information 
may  be  utilized  to  constrain  the  optimization  processes.  Three 
cases  are  included  to  illustrate  the  effects  of  contrast  media  in 
Section  V,  and  a  dozen  of  clinical  cases  without  using  contrast 
media  are  used  for  evaluation  in  Section  VI,  demonstrating  the 
effectiveness  of  our  technique.  We  conclude  in  Section  VII. 

II.  Perfusion  signals 

In  [10],  for  a  sequence  of  chest  images  I(x,y,t)  with  0  < 
/  <  255, 1  <  x  <  width  (192  for  ventilation  and  384  for 
perfusion),  1  <  y  <  height  (144  for  ventilation  and  288  for 
perfusion),  and  t  is  a  discrete  time  point  in  [0,  examtime]  (4.32 
seconds  for  ventilation  and  2.04  seconds  for  perfusion),  we  have 
defined  a  lung  functional  signal  (i.e.,  an  observation)  as  the  av¬ 
erage  pixel  intensity  of  a  region  of  interest  (ROI)  over  time: 

^2x,yeROI 

\ROI\  '  1  j 

where  \ROI\  is  the  number  of  the  pixels  in  the  region  of  inter¬ 
est.  When  the  patient  breathes  naturally,  an  observation  includes 
both  ventilation  and  perfusion  components  plus  noise.  Here,  we 
are  only  interested  in  the  perfusion  component.  Therefore,  the 
patient  is  asked  to  hold  the  breath  to  effectively  remove  the  ven¬ 
tilation  component.  For  convenience,  an  observation  in  case  of 
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Fig.  1.  A  case  with  the  breath  held  and  an  intravenous  bolus  of  X-ray  contrast 
media,  (a)  An  ROI  (region  of  interest)  in  the  right  lung  field  and  (b)  its  cor¬ 
responding  observation  -  an  enhanced  lung  perfusion  signal,  which,  due  to  the 
X-ray  physical  property,  reflects  the  blood  flow  in  the  corresponding  lung  area 
with  contrast  media. 


(a)  (b) 


Fig.  2.  A  case  with  the  breath  held  but  no  X-ray  contrast  media,  (a)  An  ROI 
(region  of  interest)  in  the  right  lung  and  (b)  its  corresponding  observation  - 
a  perfusion  signal  reflecting  the  blood  flow  in  the  lung  area  due  to  the  X-ray 
physical  property.  It  is  plotted  in  the  same  scale  as  in  Fig.  1  for  comparison. 

the  breath  held  is  called  a  perfusion  signal.  The  perfusion  sig¬ 
nal  strength  can  be  enhanced  with  an  intravenous  bolus  of  X- 
ray  contrast  media  as  shown  in  Fig.  1.  For  comparison,  Fig.  2 
shows  a  perfusion  signal  without  using  contrast  media  in  the 
same  scale. 

Pulmonary  perfusion  analysis  is  to  extract  useful  perfusion 
parameters  from  perfusion  signals  for  showing  pulmonary  per¬ 
fusion  abnormalities.  It  is  similar  to  ventilation  analysis;  conse¬ 
quently,  the  method  we  have  developed  for  ventilation  analysis 
would  be  useful  to  perfusion  analysis.  However,  perfusion  has 
different  properties  which  should  be  taken  into  consideration. 

III.  Model-based  perfusion  analysis  with  truncated 
pyramid  images 


Fig.  3.  A  perfusion  model  with  five  free  primitive  parameters:  amplitude  A 
(perfusion  strength  in  the  lung  area),  downtime  D  (time  for  the  systolic  phase  in 
the  lung  area),  uptime  U  (time  for  the  diastolic  phase  in  the  lung  area),  time  shift 
S  (time  from  the  first  image  to  the  completion  of  the  first  diastolic  phase)  and 
level  L  (the  mean  intensity  but  with  no  well-defined  medical  meaning).  This 
function  models  the  blood  change  of  pulmonary  perfusion;  it  increases  during 
the  diastolic  phase  and  decreases  during  the  systolic  phase.  The  systolic  phase 
is  generally  shorter  than  the  diastolic  phase.  The  free  parameters  downtime  D 
and  uptime  U  may  be  further  constrained  with  the  cardiac  systolic  and  diastolic 
phases  extracted  from  the  heart  to  make  the  optimization  process  fast  and  robust 
(see  Section  IV). 

as  the  ventilation  model  given  in  [10,13]  for  ventilation  analy¬ 
sis.  However,  the  five  characteristic  parameters  have  different 
medical  meanings: 

•  Amplitude  A:  Perfusion  strength. 

•  Uptime  U:  Time  corresponding  to  the  diastolic  phase  in  the 
lung  area. 

•  Downtime  D:  Time  corresponding  to  the  systolic  phase  in 
the  lung  area. 

•  Timeshift  S:  Time  from  the  first  image  to  the  completion  of 
the  first  diastolic  phase. 

•  Level  L:  Intensity  mean  -  a  mathematically  necessary  pa¬ 
rameter  without  well-defined  medical  meaning  (i.e.,  its 
value  depending  on  many  factors). 

Once  a  perfusion  model  is  available,  the  perfusion  param¬ 
eters  can  be  extracted  via  the  nonlinear  least  squares  opti¬ 
mization  procedure.  The  standard  algorithm  for  the  non¬ 
linear  least  squares  optimization  is  the  Levenberg-Marquardt 
method  [7,14,16].  This  model  is  sophisticated  enough  in  cov¬ 
erage  of  both  systolic  and  diastolic  phases  but  remains  simple 
enough  for  efficient  model  realization. 


A.  A  perfusion  model 

For  the  extraction  of  perfusion  information  from  a  perfusion 
signal,  similarly  to  ventilation  analysis,  due  to  the  small  num¬ 
ber  of  samples  and  even  poorer  signal-to-noise  ratio,  general 
signal  processing  techniques  do  not  work  well  in  the  sense  of 
result  accuracy.  Therefore,  the  model-based  approach  becomes 
more  prominent.  The  perfusion  model  consists  of  two  sinusoidal 
functions  connected  to  each  other  in  order  to  model  the  non¬ 
symmetry  of  systolic  and  diastolic  phases  as  depicted  in  Fig.  3. 
The  model  can  be  expressed  as  a  mathematical  function: 

M(A,  D,  U,  S,L,t) 

A  cosfKt'/D)  +  L  if  0  <  t'  <  D 

Acos(7r(t'  -  D)/U  +  tt)  +  L  if  D  <  t'  <  (D  +  U)  ^ 

where 

t’  =  (t-S )  mod  (D  +  U),  (3) 

and  t  indicates  time.  Please  note  that  t'  is  always  in  the  interval 
[0,D  +  U).  Mathematically,  this  perfusion  model  is  the  same 


B.  Perfusion  properties 

When  blood  flows  in  the  lungs,  the  phase  (i.e.,  timeshifts)  of  a 
pulse  signal  depends  on  its  location.  For  observations  with  dif¬ 
ferent  timeshifts,  the  average/median  operator  will  reduce  the 
amplitude  and  generate  some  artifacts.  In  other  words,  the  am¬ 
plitude  of  an  observation  corresponding  to  a  larger  lung  area  is 
smaller  than  that  of  an  observation  corresponding  to  its  subpart 
and  it  is  also  more  difficult  to  model  the  behavior  of  an  observa¬ 
tion  corresponding  to  a  larger  lung  area. 

For  simplicity,  let’s  assume  that  two  adjacent  subareas  have 
two  pulse  signals  which  are  perfect  sinusoidal  functions  with 
the  same  amplitude  and  frequency  but  different  initial  phases: 

Pi  =  A  sin(27r/£  +  Of)  (4) 

P2  =  A  sin(27 rft  +  02)  (5) 

where  A  is  the  perfusion  amplitude,  /  is  the  pulse  frequency, 
and  Oi  and  62  are  the  initial  phases  of  the  two  pulse  signals,  re¬ 
spectively.  According  to  the  basic  properties  of  trigonometrical 
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sum,  the  average  pulse  for  the  total  area  is 

P  =  Oi+P2)/2  (6) 

=  A  sin[27r/£  +  (0i  +  02)/2]  cos[(0i  -  02)/2]  (7) 

=  {A  cos[(0i  -  02 ) / 2] }  sin{27r/t  +  [(0i  +  02)/2]}(8) 

This  naturally  implies  that  as  long  as  (0i  —  02)/2  ^  2tt  •  iV, 
where  N  =  ***,  — 1,0,1,***,  the  amplitude  of  the  pulse  signal  p 
(a  new  observation)  will  be  smaller  than  that  of  the  two  original 
pulse  signals  p\  and  p2,  and  its  phase  (the  average  of  0i  and  02) 
also  differs  from  those  of  pi  and  p2 . 

If  the  pulse  signals  were  perfect  sinusoidal  functions  with  the 
same  frequency,  we  would  have  another  perfect  sinusoidal  func¬ 
tion  with  the  same  frequency  by  the  average  operator  as  shown 
above,  although  it  may  have  a  smaller  amplitude  and  a  different 
initial  phase.  We  may  reasonably  assume  that  all  pulse  signals  of 
a  patient  during  the  short  examination  have  the  same  frequency, 
but  the  reality  is  that  the  pulse  signals  are  not  perfect  sinusoidal 
functions  (/.<?.,  the  systolic  and  diastolic  phases  are  not  symmet¬ 
rical).  Therefore,  the  average  operator  not  only  reduces  the  am¬ 
plitude  but  also  generates  some  artifacts,  since  we  can  not  guar¬ 
antee  that  a  set  of  the  parameters  (C",  A' ,  D'  ,U' ,  S' ,  L')  satisfy 
the  following  equation: 

M(A,  D,  U,  Si,  L,  t)  +  M(A,  D,  U,  S2,  L,  t)  (9) 
=  C'M(A,,D,,U',S,1L,A) 

provided  that  D  ^  U  and  the  difference  of  Si  and  S2  is  not  a 
cycle  length  (i.e.,  U  +  D).  In  other  words,  this  makes  it  diffi¬ 
cult  to  model  the  behavior  of  an  observation  corresponding  to 
a  larger  lung  area.  We  have  demonstrated  this  idea  in  [8]  with 
artificial  signals:  The  average/median  operator  reduces  the  re¬ 
sultant  signal’s  amplitude  and  also  generates  some  artifacts  {i.e., 
the  model  cannot  be  exactly  fitted  to  the  resultant  signals  from 
the  average  or  median  operator). 

This  means  that  for  perfusion  analysis,  we  need  to  use  higher 
spatial  resolution  images  in  order  to  avoid  the  amplitude  reduc¬ 
tion  and  the  artifacts. 

C.  Truncated  pyramids 

Just  like  ventilation  analysis  [10,13],  we  are  interested  in  not 
only  the  overall  perfusion  function  of  a  lung  segment  but  also  its 
subsegments.  This  naturally  leads  to  multiresolutional  analysis 
based  upon  pyramids.  The  pyramid  for  each  image  is  obtained 
by  recursively  dividing  the  lung  field  into  4  (approximately) 
equal  regions  till  the  resulting  region  is  too  small  to  be  divided; 
the  intensity  value  of  a  pixel  (region)  in  the  pyramid  is  calculated 
as  the  average  intensity  of  its  corresponding  area.  This  process 
results  in  a  pyramid  of  observations.  However,  because  of  the 
perfusion  properties  discussed  in  Section  III-B,  the  higher  lev¬ 
els  of  the  pyramids  should  be  ignored  due  to  their  poorer  spatial 
resolutions.  This  yields  a  truncated  pyramid  where  the  lowest 
spatial  resolution  is  still  sufficiently  high  to  avoid  the  amplitude 
reduction  and  the  artifacts.  The  truncated  pyramid  is  processed 
level  by  level  in  a  top-down  fashion  to  extract  the  perfusion  in¬ 
formation  from  each  of  its  observations. 


IV.  Robustly  accelerating  perfusion  analysis  with  cardiac 
information 

We  have  formulated  the  extraction  of  perfusion  parameters 
as  a  nonlinear  least  squares  optimization  problem.  It  is  well 
known  that  the  convergence  speed  and  result  accuracy  depend 
on  the  initial  guess.  Therefore,  it  is  essential  to  have  a  good 
guess  when  fitting  the  model  to  an  observation.  The  truncated 
pyramid  can  partially  help  in  this  end,  since  the  optimized  re¬ 
sult  of  an  observation  in  the  truncated  pyramid  can  be  taken  as 
an  initial  guess  for  its  immediate  child  observations.  However, 
the  top  level  needs  to  be  handled  separately.  Due  to  its  small 
signal-to-noise  ratio,  it  is  difficult  to  estimate  an  initial  guess  for 
a  perfusion  signal.  It  is  also  the  small  signal-to-noise  ratio  that 
gives  more  local  minima  for  the  error  function  in  optimization, 
consequently,  it  takes  more  time  to  converge  to  a  solution.  Nat¬ 
urally,  it  would  be  desirable  if  we  can  reduce  the  number  of  free 
perfusion  model  parameters,  since  it  not  only  reduces  the  opti¬ 
mization  time  but  also  improves  its  stability  and  result  accuracy. 
In  the  following,  we  shall  present  a  simple  yet  accurate  approach 
to  extract  cardiac  systolic  and  diastolic  phases  from  the  heart,  so 
that  this  cardiac  information  may  be  utilized  to  constrain  the  op¬ 
timization  process,  making  perfusion  analysis  not  only  fast  but 
also  robust. 

A.  Extracting  systolic  and  diastolic  phases  from  the  heart 

The  perfusion  examination  takes  only  2-3  seconds,  thus,  it 
would  be  reasonable  to  assume  that  the  duration  of  the  patient’s 
systolic  phase  is  the  same  anywhere  in  the  lungs  during  the  ex¬ 
amination  and  so  is  the  duration  of  the  diastolic  phase.  Conse¬ 
quently,  the  patient’s  pulse  frequency  is  the  same  anywhere  in 
the  lungs  during  the  examination.  Based  on  this  assumption,  we 
can  extract  the  systolic  and  diastolic  phases  from  one  source: 
the  heart,  so  that  the  estimated  results  of  the  systolic  and  dias¬ 
tolic  phases  can  be  used  as  fixed  parameters  to  accelerate  the 
convergence  of  the  optimization  processes. 

More  specifically,  first  we  employ  a  trick  by  using  an  ROI 
on  the  heart  border  as  shown  in  Fig.  4(a)  to  have  an  observa¬ 
tion  (also  called  a  heart  signal)  (see  Fig.  4(b)).  The  dominant 
information  this  observation  carries  is  the  change  of  the  heart 
proportion  in  the  ROI  from  one  frame  to  another.  This  signal  is 
generally  strong,  so  the  initial  guesses  for  those  parameters  can 
be  conveniently  estimated  from  the  signal  itself.  The  uptime  of 
this  signal  corresponds  to  the  systolic  phase  of  the  heart,  while 
its  downtime  corresponds  to  the  diastolic  phase  of  the  heart.  By 
fitting  the  perfusion  model  to  this  observation,  the  systolic  and 
diastolic  phases  are  available.  Mathematically,  from  the  heart 
observation  Oh(t),  the  fitting  can  determine  a  set  of  parameters 
{A*h,  and  L£)  which  minimize 

Y  [M{Ah,Dh,Uh,Sh,Lh,t)  -  Oh{t)f.  (10) 

tE[0,examtime\ 

With  the  United  Snakes  technique  [12,1 1],  we  have  justified  that 
this  simple  trick  is  actually  accurate  enough  in  extracting  sys¬ 
tolic  and  diastolic  phases  for  pulmonary  perfusion  analysis  [8]. 
Alternatively,  we  may  use  an  ROI  within  the  heart  area  to  ex¬ 
tract  the  systolic  and  diastolic  phases.  The  problem  is  that  such 
an  observation  is  rather  noisy  and  weak.  Moreover,  it  is  com- 
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(a) 


(b) 


Fig.  4.  An  example  for  extracting  the  systolic  and  diastolic  phases  from  the 
heart  on  the  case  with  contrast  media,  (a)  An  ROI  (region  of  interest)  on  the 
heart  border,  (b)  The  corresponding  observation  and  the  parameter  extraction 
process.  The  observation  indicated  with  “o”  mainly  reflects  the  change  of  the 
heart  proportion  in  the  ROI  from  frame  to  frame.  The  initial  guess  is  plotted  as 
dashed  curve  and  the  final  solution  is  the  solid  curve.  During  the  systolic  phase, 
the  heart  proportion  in  the  ROI  becomes  smaller  and  smaller,  thus,  the  average 
intensity  values  of  the  ROI  gets  bigger  and  bigger.  In  other  words,  the  uptime 
of  this  signal  corresponds  to  the  systolic  phase  of  the  heart;  while  its  downtime 
corresponds  to  the  diastolic  phase  of  the  heart  —  the  uptime  and  downtime  ex¬ 
tracted  from  a  heart  signal  have  completely  different  medical  meanings  from 
those  of  an  observation  in  the  lung  (see  Fig.  5).  The  medical  meaning  of  the 
extracted  amplitude  from  the  heart  signal  is  undefined,  since  not  only  does  it 
depend  on  the  heart  pumping  strength  but  also  the  amount  of  air  in  the  lungs. 

plicated  in  its  medical  meaning  because  of  the  overlapping  lung 
area. 

B.  Constraining  the  fitting  process 

The  estimated  systolic  and  diastolic  phases  (U%  and  D f)  from 
the  heart  signal  are  then  used  as  fixed  parameters  in  extracting 
the  perfusion  parameters  from  observations  in  the  lung  fields 
(see  Fig.  5).  However,  it  should  be  noted  that  the  uptime  and 
downtime  of  an  observation  in  the  lung  fields  have  completely 
different  meanings  from  those  of  the  heart  signal:  the  downtime 
of  the  signal  in  the  lung  areas  corresponds  to  the  systolic  phase; 
while  its  uptime  corresponds  to  the  diastolic  phase.  Mathemati¬ 
cally,  from  the  lung  signal  0/  (t) ,  we  determine  a  set  of  parame¬ 
ters  (Ai ,  Di  ,  Ui  ,  Si  and  L z*)  which  minimize 

Y,  [M{AuDi,UuSuLi,t)-Oi{t)f  (11) 

t£[0,examtime\ 

subject  to  the  constraints: 

Di  =  U*hJ  Ut  =  D*h.  (12) 

Naturally,  we  always  have 

d;  =  ui  u;  =  D*h.  d3) 

Therefore,  the  perfusion  analysis  gives  three  parameter  images 
(perfusion  amplitude  image,  timeshift  image  and  level  image). 
As  mentioned  before,  the  parameter  L  (level)  has  no  well- 
defined  medical  meaning.  Furthermore,  we  could  not  see  any 
clear  perfusion  waves  in  the  lung  fields  from  the  timeshift  image. 
We  speculate  that  this  is  due  to  the  low  spatial  resolution  used 
in  the  current  research.  Therefore,  we  only  use  the  perfusion 
amplitude  image  which  appears  to  be  sufficient  for  diagnosis  of 
pulmonary  embolism. 

V.  Effects  of  contrast  media 

Three  clinical  cases  have  been  used  as  experiments  for  per¬ 
fusion  analysis  with  contrast  media.  Due  to  historical  reasons, 


Fig.  5.  Using  the  cardiac  systolic  and  diastolic  phases  to  constrain  the  parameter 
extraction  from  an  enhanced  pulmonary  perfusion  signal  with  contrast  media, 
(a)  An  ROI  in  the  right  lung  field,  (b)  The  perfusion  signal  (indicated  with  “o”) 
and  the  parameter  extraction  process  (the  dashed  curve  for  the  initial  guess  and 
the  solid  curve  is  the  final  solution).  The  downtime  of  a  pulmonary  perfusion 
signal  corresponds  to  its  systolic  phase  (more  blood  in  the  lung  area);  while  the 
uptime  corresponds  to  its  diastolic  phase  (less  blood  in  the  lung  area). 

the  image  matrix  is  192  x  144,  where  a  pixel  corresponds  to  an 
area  of  roughly  3  mm  x  3  mm  of  the  lung.  It  is  the  same  spatial 
resolution  as  we  have  used  in  ventilation  analysis  [10,13],  but 
the  temporal  sampling  frequency  has  been  doubled  (i.e.,  25  Hz). 

Based  on  the  three  perfusion  amplitude  images  in  Fig.  6,  con¬ 
trast  media  concentration  in  pulmonary  vessels  appears  to  have 
been  constant  in  Case  1 ,  while  in  Case  2  there  is  inflow  of  con¬ 
trast  media  into  pulmonary  arteries  causing  strong  arterial  sig¬ 
nal  and  heart  motion  has  caused  broader  motion  artifact  than  in 
Case  1 .  Case  3  represents  the  inflow  period  of  contrast  media 
through  the  right  subclavian  vein  and  pulmonary  arteries  with 
little  signal  change  in  pulmonary  veins  ( e.g .,  early  phase  after 
contrast  media  injection). 

This  experiment  has  demonstrated  that  contrast  media  can 
significantly  improve  the  pulmonary  perfusion  signal  strength, 
but  may  cause  some  artifacts  disturbing  the  parameter  image  in¬ 
terpretation.  Furthermore,  contrast  media  are  expensive,  carry  a 
risk  of  contrast  media  reactions,  should  not  be  used  in  patients 
with  pulmonary  edema  or  any  renal  problem,  and  also  require 
timing  in  taking  the  X-ray  series.  Therefore,  no  contrast  media 
will  be  used  in  clinical  evaluation. 

VI.  Clinical  case  studies 

In  clinical  evaluation,  we  have  selected  18  patients,  who  were 
referred  to  this  examination  mainly  to  exclude  pulmonary  em¬ 
bolism.  All  of  them  were  examined  with  no  contrast  media.  The 
image  resolution  is  384  x  288  pixels  and  the  sampling  frequency 
is  25  Hz.  In  the  following,  we  shall  classify  the  perfusion  ab¬ 
normal  findings  into  three  types  and  present  three  representative 
cases. 

A.  Three  types  of  perfusion  abnormalities 

We  have  classified  perfusion  abnormal  findings  into  the  fol¬ 
lowing  three  types: 

•  No  perfusion  (NP):  The  perfusion  amplitude  is  extremely 
small  or  zero.  This  is  associated  with  the  complete  occlu¬ 
sion  of  pulmonary  arteries  by  an  embolism,  and  often  seen 
in  the  upper  and  middle  parts  of  the  lung.  When  the  no¬ 
perfusion  area  becomes  larger,  likely  there  will  be  associ¬ 
ated  overactive  perfusion  (OP)  in  the  normal  part(s)  of  the 
lungs. 

•  Reduced  perfusion  (RP):  The  perfusion  amplitude  becomes 
smaller  than  expected  and  can  be  seen  as  darker  areas  in 
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(a)  (b)  (c) 

Fig.  6.  Perfusion  amplitude  images  of  three  clinical  cases  with  contrast  media:  (a)  contrast  media  concentration  in  this  case  is  uniformly  distributed,  (b)  there  is 
an  inflow  of  contrast  media  into  pulmonary  arteries  causing  strong  arterial  signal  (indicated  with  an  arrow),  (c)  This  case  represents  the  inflow  period  of  contrast 
media  through  the  right  subclavian  vein  and  pulmonary  arteries  indicated  with  arrows. 


(a)  Perfusion  amplitude  (b)  CT  slice  1  (c)  CT  slice  2 


Fig.  7.  Case  1:  (a)  Overactive  perfusion  (OP)  is  seen  in  the  right  lower  lung  field  (indicated  with  a  bracket),  while  slightly  reduced  perfusion  (RP)  is  shown  in  the 
rest  of  the  right  lung.  Reduced  perfusion  (RP)  is  revealed  in  the  central  part  of  the  left  lung  (indicated  with  an  arrow).  The  CT  images  (b,  c)  of  the  same  patient 
show  embolic  masses  partially  filling  the  right  pulmonary  artery  and  also  material  in  the  left  lower  lobe  artery  (indicated  with  arrows).  Our  findings  also  show  a 
good  correlation  with  the  patient’s  clinical  report  of  the  NM  study  performed  in  another  hospital  (the  NM  image  is  not  available). 


the  perfusion  amplitude  image.  This  is  the  typical  phe¬ 
nomenon  of  pulmonary  embolism  with  partial  occlusion. 

•  Overactive  perfusion  (OP):  The  perfusion  amplitude  is  big¬ 
ger  than  expected  and  the  area  should  be  considered  as  nor¬ 
mal.  This  is  the  phenomenon  caused  by  the  excessive  blood 
flow  redirected  into  the  normal  area  due  to  no-perfusion 
and  reduced  perfusion  in  other  parts  of  the  lungs. 

The  lung  field  may  be  visually  divided  into  three  regions  (up¬ 
per,  middle  and  lower).  The  statistics  of  perfusion  abnormalities 
found  in  the  1 8  patients  are  summarized  in  Table  1 . 


Right/Left 

NP 

RP 

OP 

Upper 

2/3 

5/5 

1/1 

Middle 

2/1 

3/3 

0/9 

Lower 

1/2 

7/2 

1/3 

Table  1.  Statistics  of  perfusion  abnormalities. 

B.  Three  representative  clinical  cases 

Illustrated  here  are  three  representative  cases,  who  were  also 
examined  with  computed  tomography  (CT)  and  pulmonary  per¬ 
fusion  nuclear  medicine  (NM).  Both  CT  and  NM  are  the  “golden 
standard”  method  in  detection  of  pulmonary  perfusion  distur¬ 
bances  ( e.g .,  pulmonary  embolism).  NM  shows  the  distribu¬ 
tion  of  pulmonary  perfusion,  while  CT  reveals  the  thrombotic 
masses  causing  pulmonary  embolism. 

•  Case  1  (see  Fig.  7):  Pulmonary  embolism  of  the  right  mid¬ 
dle  lobe  and  the  right  upper  lobe  is  associated  with  reduced 
perfusion  (RP)  in  the  middle  and  upper  fields  of  the  right 
lung.  In  addition,  there  is  reduced  perfusion  (RP)  in  the  left 
upper  lobe  and  perihilar  region  of  the  left  lower  lobe.  The 
reason  for  a  very  high  amplitude  (overactive  perfusion,  OP) 
in  the  right  lower  lung  field  is  due  to  the  high  concentration 
of  the  blood  in  this  area.  These  findings  show  a  good  cor¬ 
relation  with  both  CT  and  NM  studies. 

•  Case  2  (see  Fig.  8):  Pulmonary  embolism  in  the  right  lung 


(a)  Perfusion  amplitude  (b)  Nuclear  medicine 


Fig.  8.  Case  2:  (a)  Overall  reduced  perfusion  (RP)  of  the  right  lung.  Overactive 
perfusion  (OP)  is  seen  centrally  in  the  left  lung,  while  reduced  perfusion  (RP)  is 
shown  in  the  left  apex  and  a  small  area  (indicated  with  an  arrow)  in  the  upper  left 
lung  field.  Pulmonary  embolism  in  the  right  lung  and  in  the  superior  segment 
of  the  left  lower  lobe  shown  by  CT  and  NM  (b)  studies.  Generally,  the  NM 
image  shows  the  perfusion  activity  in  the  anterior  parts  of  the  lungs,  while  our 
perfusion  amplitude  image  reveals  perfusion  through  the  lungs.  Therefore,  the 
reduced  perfusion  of  the  lower  part  of  the  right  lung  shown  by  our  method  (a) 
may  be  explained  by  the  thrombotic  mass  in  the  hilar  region  of  the  right  lung  as 
reported  by  CT. 

and  in  the  superior  segment  of  the  left  lower  lobe  shown  by 
CT  and  NM  studies  is  correlated  to  reduced  perfusion  (RP) 
of  the  right  lung  and  reduced  perfusion  (RP)  of  the  left  apex 
and  a  small  area  in  the  upper  left  lung  field.  Overactive 
perfusion  (OP)  is  seen  centrally  in  the  left  lung. 

•  Case  3  (see  Fig.  9):  Generally  reduced  perfusion  (RP)  of 
the  right  lung  and  no  perfusion  (NP)  of  the  lateral  recess 
in  the  left  lung.  These  findings  are  consistent  with  CT  and 
NM  studies.  Overactive  perfusion  (OP)  seen  centrally  in 
the  left  lung  is  due  to  the  redirection  of  blood  flow. 

The  clinical  case  studies  show  that  our  model-based  method 
for  pulmonary  perfusion  analysis  is  effective  for  pulmonary  em¬ 
bolism  even  without  using  contrast  media,  demonstrating  con¬ 
sistent  correlations  with  CT  and  NM  studies.  This  gives  our 
present  technique  some  advantages.  It  takes  only  about  2  sec¬ 
onds  and  involves  no  radioactive  isotopes,  no  contrast  media  and 
only  low  radiation  dose  [6,5,17].  The  NM  study  takes  much 
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(a)  Perfusion  amplitude  (b)  Nuclear  medicine 


Fig.  9.  Case  3:  (a)  Generally  reduced  perfusion  (RP)  of  the  right  lung.  No 
perfusion  (NP)  of  the  lateral  recess  in  the  left  lung  (indicated  with  an  arrow). 
Overactive  perfusion  (OP)  seen  centrally  in  the  left  lung.  These  findings  are 
consistent  with  CT  and  NM  (b)  studies. 


(a)  (b) 

Fig.  10.  Case  3:  (a)  Perfusion  amplitude  image  and  (b)  perfusion  frequency 
image  obtained  from  an  analysis  without  using  cardiac  information.  Compared 
with  the  amplitude  image  in  Fig.  9(a),  this  amplitude  image  (a)  is  rather  noisy  in 
the  no-perfusion  and  reduced  perfusion  areas.  The  no-perfusion  area  of  the  lat¬ 
eral  recess  in  the  left  lung  is  not  so  visible  as  that  in  Fig.  9(a).  The  no-perfusion 
and  reduced  perfusion  areas  in  the  perfusion  frequency  image  (b)  are  also  noisy. 

longer  (over  20  minutes)  and  it  is  not  readily  available  in  any 
hospital.  The  CT  study  uses  high  amount  of  contrast  media  and 
has  side  effects.  Furthermore,  CT  is  expensive  and  cannot  be 
used  for  all  patients  (e.g.,  anxiety  and  contrast  media  reactions). 

This  performance  of  our  computer  analysis  method  is  the  con¬ 
sequence  of  the  idea  of  reducing  the  number  of  free  perfusion 
model  parameters  (see  Section  IV),  which  makes  the  optimiza¬ 
tion  process  not  only  fast  but  also  robust  in  computation.  For 
comparison,  we  have  performed  a  test  on  Case  3  without  con¬ 
straining  uptime  and  downtime  with  the  cardiac  information.  It 
takes  about  10  times  longer  and  the  results  are  rather  sensitive  to 
noise,  specially  in  the  no-perfusion  and  reduced  perfusion  areas 
as  shown  in  Fig.  10.  Perfusion  analysis  is  to  show  pulmonary 
perfusion  abnormalities  (i.e.,  no-perfusion  and  reduced  perfu¬ 
sion).  In  a  reduced  perfusion  area,  the  signal  component  with  the 
pulse  frequency  is  weak,  while,  in  a  no-perfusion  area,  its  obser¬ 
vation  has  no  such  a  component.  When  the  uptime  and  down¬ 
time  are  constrained,  the  extraction  process  will  only  search  for 
the  component  with  the  pulse  frequency  in  the  observation,  con¬ 
sequently,  it  will  be  able  to  quickly  and  robustly  give  a  small 
value  or  zero  to  the  amplitude  parameter  in  the  reduced  perfu¬ 
sion  area  or  no-perfusion  area  without  being  disturbed  by  noise. 

VII.  Conclusion 

We  have  extended  our  method  for  ventilation  analysis  to  pul¬ 
monary  perfusion  analysis.  According  to  perfusion  properties, 
we  have  assigned  new  medical  meanings  to  our  model  parame¬ 
ters  in  order  to  form  a  perfusion  model  and  truncated  the  pyra¬ 
mid  structure  so  as  to  avoid  perfusion  artifacts. 

Three  clinical  cases  have  been  used  to  illustrate  the  effects  of 
contrast  media.  To  enable  comparison  with  computed  tomogra¬ 
phy  (CT)  and  nuclear  medicine  (NM)  studies,  perfusion  abnor¬ 


mal  findings  are  classified  into  three  types.  The  clinical  eval¬ 
uation  has  shown  that  our  computer  analysis  method  is  effec¬ 
tive  for  pulmonary  embolism  even  without  using  contrast  media, 
demonstrating  consistent  correlations  with  CT  and  NM  studies. 
This  performance  is  the  consequence  of  the  idea  that  the  cardiac 
information  recorded  in  the  perfusion  image  sequence  may  be 
utilized  to  accelerate  pulmonary  perfusion  analysis  and  improve 
its  sensitivity  in  detecting  pulmonary  embolism.  In  doing  so,  a 
simple  yet  accurate  approach  has  been  introduced  to  extract  car¬ 
diac  systolic  and  diastolic  phases  from  the  heart  for  constrain¬ 
ing  the  optimization  processes.  This  idea  has  not  only  made 
perfusion  analysis  fast  but  also  robust;  consequently,  perfusion 
analysis  becomes  feasible  without  using  contrast  media. 
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